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We present an analytical study of three-state Feshbach resonances induced by second-order cou- 
plings. Such resonances arise when the scattering amplitude is modified by the interaction with a 
bound state that is not directly coupled to the scattering state containing incoming flux. Coupling 
occurs indirectly through an intermediate state. We consider two problems: (i) the intermediate 
state is a scattering state in a distinct open channel; (ii) the intermediate state is an off-resonant 
bound state in a distinct closed channel. The first problem is a model of electric-field-induced reso- 
nances in ultracold collisions of alkali metal atoms [Phys. Rev. A 75, 032709 (2007)] and the second 
problem is relevant for ultracold collisions of complex polyatomic molecules, chemical reaction dy- 
namics, photoassociation of ultracold atoms, and electron - molecule scattering. Our analysis yields 
general expressions for the energy dependence of the T-matrix elements modified by three-state res- 
onances and the dependence of the resonance positions and widths on coupling amplitudes for the 
weak-coupling limit. We show that the second problem can be generalized to describe resonances in- 
duced by indirect coupling through an arbitrary number of sequentially coupled off-resonant bound 
states and analyze the dependence of the resonance width on the number of the intermediate states. 



I. INTRODUCTION 



n 

Following the paper of Tiesinga, Verhaar and Stoof [l| , Feshbach resonances have been used as an important tool 
for controlling interactions of ultracold atoms 2fl , the creation of ultracold molecules 3|, U, la, Ifj, LZ| and experimental 



studies of correlated phenomena in ultracold gases 8|, [9 



101 ] . A Feshbach resonance occurs when a scattering state of 



two colliding particles interacts with a metastable bound state of the two-particle system. The coupling between the 
scattering state and the bound state leads to a resonant enhancement of the scattering cross section as the energy 
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of the scattering state approaches the energy of the bound state. Collisions at ultracold temperatures are entirely 
determined by single partial-wave scattering: s-wave scattering for collisions of bosons or distinct atoms and p-wave 
scattering for collisions of identical fermions. Feshbach resonances in ultracold collisions can therefore be described 



by 



a two-state model where a single partial-wave state interacts with an isolated bound state ( 3j, 



ll| and references 



llj). The energy of the bound state can be tuned with respect to the collision threshold by a magnetic field so 



3 



magnetic fields have been used in most experiments on Feshbach resonances in ultracold gases ([12| and references 
therein). We have recently proposed that Feshbach resonances in ultracold collisions of distinct atoms may also be 



tuned by dc electric fields [l3j. An electric field induces couplings between different partial waves of the colliding atoms, 
which results in a different mechanism of Feshbach resonances. For example, it is not necessary that the dominant 
scattering state be directly coupled to the bound state giving rise to resonances so the resonant enhancement of the 
scattering cross-section may be determined by second order couplings and the resonances cannot be described by a 
two-state model. The purpose of this paper is to develop a formal theory for three-state Feshbach resonances induced 
by second-order couplings. 



14 



la ] to analyze three-state Feshbach 



In the first part of this paper, we use the projection-operator method 
resonances induced by indirect couplings via an intermediate scattering state. Figure [1^ illustrates the mechanism of 
such resonances. The initial scattering state a is coupled to another scattering state (3, while [3 is coupled to bound 
state 7. The sequence of a-f3 and [3-j couplings gives rise to a resonant enhancement in the state a. This is a model 



for electric-field-induced resonances in ultracold collisions of distinct atoms, such as Li and Cs, described in Refs 
and 



al 



13] 



17j . Diirr et 



16j. Three-state resonances of this type can also be induced by magnetic dipole-dipole interaction 
171 ] have recently presented a study of magnetic-dipole-induced predissociation of ultracold Rb2 molecules leading 
to population of multiple scattering states. 

Resonances induced by indirect couplings are also important for collisions of large polyatomic molecules with 
ultracold atoms. One proposed method of cooling polyatomic molecules to ultracold temperatures is sympathetic 
cooling. In the absence of reaction processes, molecules can potentially be cooled by elastic collisions in a reservoir of 



ultracold atoms 



18 



19 



20( . The experimental realization of this method may, however, be complicated by naturally 
occurring Feshbach resonances. The density of molecular states near zero point energy is very large in polyatomic 
molecules [21| and these bound states may give rise to Feshbach resonances in collisions between molecules and 
ultracold atoms. The probability of three-body recombination and other loss processes in collisions involving large 
polyatomic molecules may thus be enhanced. Not all of the bound states may, however, give rise to Feshbach 



3 
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FIG. 1: Schematic diagram of a three-state Feshbach resonance in channel a induced by indirect couplings to a bound state in 
channel 7 via channel /3: in (a) /3 is a scattering state in an open channel; (b) j3 is a bound state in a closed channel. 

resonances. The interaction between an atom and a large molecule usually probes only a limited number of molecular 
states. The other molecular states can interact with the atom - molecule scattering state by indirect couplings and it 
is not clear whether they can generate Feshbach resonances with significant widths. To elucidate the role of molecular 
states indirectly coupled to the scattering state, we analyze Feshbach resonances mediated by a sequence of one 
continuum-bound and several bound-bound couplings. We first consider a three-state Feshbach resonance induced 
by a combination of two couplings: a coupling between the scattering state and an off-resonant bound state and a 
coupling between the off-resonant bound state and a resonant bound state (Fig.[lJ>). In the third part of the paper, 
we generalize this analysis to a system with Feshbach resonances induced by higher-order couplings. 

Bohn and Julienne examined similar Feshbach resonances, with second- and higher-order couplings through a series 



22j . Their derivation yields the 



of bound states, in collisions of ultracold atoms induced by optical laser fields 
dependence of the scattering amplitude on the Rabi frequency and the detuning of the laser fields. The problems 
considered in Sections IIVI and M of this paper differ from their systems in that, for cold atom-molecule collisions, the 
couplings between states depend on the atom-molecule distance and vanish asymptotically. 

The problems we solve in this paper are equivalent to the diagonalization of a Hamiltonian containing one or more 
bound states coupled to one or more distinct continua. Fano presented solutions for the diagonalization of a number 

n 

of similar Hamiltonians with coupled bound states and continua [231 ] . Our systems differ from the problems solved by 
Fano in the patterns of couplings present. For example, Fano considered one bound state coupled directly to multiple 
continua not coupled to each other, or a single continuum state coupled to multiple bound states with no couplings 
between bound states. Feshbach described a general method for the solution of a continuum coupled to multiple 
bound states with couplings between states ([24j|. p. 157); Sections IIVI and IVl represent an application of an approach 
similar to Feshbach's to the special case of systems of tridiagonal form. 



II. PROJECTION OPERATOR METHODS FOR FESHBACH RESONANCES 



Our approach is based on the projection operator method of Feshbach 14|, [l5j, |24| . In this section, we recapitulate 
the method in order to define the notation. The Hilbert space of the total Hamiltonian H for the collision system is 
partitioned into two orthogonal Hilbert spaces: P, comprising all open channels of the system; and Q, containing all 
closed channels. These will be referred to as the P-channel and the Q-channel, although in practice they may contain 
multiple channels. 

The total wave function of the system, \^> + (E)), belongs to the entire Hilbert space, and satisfies the Schrodinger 
equation 



(E-H)\y + (E)} = 0, 



(1) 



where E is the total energy of the colliding particles. The projections of |^ + ) onto the P and Q subspaces are given 
by 



P|*+) = |*+) 

Q|* + )^|* Q ), 



(2) 
(3) 



where P and Q are projection operators. The projected state vectors satisfy a system of coupled equations 

{E-H PP )\^+) = H PQ \* Q ) (4) 
(E - H QQ )\* Q ) = H QP \*+) , (5) 

where H PP = PHP, H QQ = QHQ, H PQ = PHQ, and H QP = QHP. 
Formally inverting Eq. (|5|) with the Green's operator (E — Hqq)~ 1 gives 

l*Q> = rp \ H QP \*+) , (6) 

which, when substituted into Eq. yields 

(E - H PP )\*+) = H PQ e ^ Hqq H qp \*+) . (7) 

This is an effective Schrodinger equation for \^p) 

(E-H eS (E))\HF+)=0, (8) 

with an energy-dependent pseudo-Hamiltonian H c g(E) = H PP + H P q(E — Hqq)~ 1 Hq P . The effective potential 
acting in the P-channel, after the effects of the coupling to the Q-channel are taken into account, is thus 

V cS (E) = V PP + H PQ 1 H QP , (9) 
h _ H QQ 

where V pp is the potential energy term of the Hamiltonian H pp . This procedure eliminates \^Sq) from the problem. 
A scattering state solution for the isolated P-channel \(pp(E)) satisfies 

(E - Hp P )\4>p) = 0, (fO) 

and a free scattering state in the P subspace \x P (E)) can be obtained from 

(E-Kpp)\xp) = 0, (11) 

where K PP is the kinetic energy operator of the Hamiltonian H PP such that H PP — K PP + V pp . The functions \xp) 
and \4>p) are related by 

1#) = \Xp) + E+ l Hpp V Pr\*r) ■ ( 12 ) 

The + superscript on the Green's operator E +^ Hpp indicates that the calculation is to be performed using the 
operator E+ill \ Hpp , with r/ a small positive real number, and the limit of i] — > + should be taken at the end of 



the calculation. Eq. (fT2|) gives the scattering state solution for an outgoing wave. We assume that the system is 
invariant under time reversal, so the incoming wave scattering state solution corresponding to a free state \xp') is 
0p/) = \xp') + (E~ — Hpp)~ 1 Vpp\xp')- The minus superscript indicates the operator E _ iT] > _ Hpp with the limit 
77 — > + taken at the end of the calculation. The subscripts P, P' on state vectors denote particular channels within 
the P subspace. 

Applying the Green's operator (E + — Hpp)^ 1 to Eq. ([7]), we obtain an implicit equation for \^f P ), 

|<E+) = |0+)+ 1 H PQ 1 H QP \#+). (13) 

— tlpp hj — tlQQ 

With Eqs. (jT2J) and (fT3"|) , our choice of the particular free state \xp) selects the incident state and the particular 
solutions \4>~p) and |^p). The term \4>~p) appears on the right-hand side of Eq. (fT3} because the general solution of 
Eq. |(7J| contains a term in the null space of the operator E — Hpp which determines the particular solution. 

We are interested in collision- induced transitions to the free state \xp>(E)), which in general may not be the same 
as \xp{E)). The adjoint of the expression for \<f> P ,) is 

(4>p>\ = (xp-\ + {xp'\Vpp e+ ^ Rpp ■ (14) 

The on-shell (i.e. between states of the same energy) T- matrix element Tprp(E) is the probability amplitude for a 
transition from the incoming free state \xp(E)} to the outgoing free state \xp>(E)}. It is given by 

Tp< P = (xp>\V eS \y+) 

= (xp'\Vpp\0+) + (d>p,\H PQ - 1 H QP \^+) (15) 

& - tlQQ 

rpb 1 rpr 

— 1 P'P T 1 P'P ) 

where T P , P is the transition amplitude for background scattering in the isolated P-channel and T P , P is the resonant 



contribution due to coupling to the Q-channel 25]. Eqs. 113} and (jT4j) have been used to obtain the second line of 
Eq. |T5}. We assume that the Q-channel contains an isolated bound state \<f>\ ) giving rise to the resonance, so that 

1 „ \<h){<h\ 



E - Hqq E — ei 

where ei is the energy of the bound state. Using this expression, we can write for T P , P , 



(16) 



rpr _ (<j> P ,\H P Q\cj) 1 ){t 1 \HQp\*+) 

T PIP - • (17) 



7 

To evaluate the matrix element (^>i\Hqp\^>^) we multipy Eq. (|T3|) on the left by (4>i\Hqp and solve to obtain 

(^\H QP \n) = < ^ l|gqp| ^ ) . (is) 

(0i I H QP e+ } Hpp H PQ \(j)x) 
1 E-e 1 

Substituting this into Eq. (fl~7|) we obtain 

T r _ (<t>p>\Hp Q \<f>Q)((f>Q\HQp\(f) p ') 

E-e x -{4> Q \H QP 

E+-H P p H Pq\4>q) 

The Green's operator in the denominator of this equation may be expanded in the basis \<j)p~„(E')), 



poo 

/ dE '^ 

Jo P „ 



E+-H P p J $t E+-E' 

which gives 



(20) 



Re ((f>Q\H Q p E+ _ Hpp H PQ\<i ) Q) =V J o dE E-E> = ' ' 



(21) 



where V indicates the Cauchy principal value integral, in which the symmetric limit approaching the singularity at 
E' = E is taken. Note that, as a result of taking the real part of the expression, the complex- valued E + in the energy 
denominator of Eq. (|20|) is replaced by the real-valued E. 
Since 

lim = 2mS(E - E') , (22) 

^o+ E-irj-E' E + irj-E' v ; ' y ' 



we obtain 



Im (4>q\Hqp * Hp Q \^ q ) = -nY,m\H QP \rt„(E))\ 2 = ™ ■ (23) 



With the additional definition 



E+-H P p 



Tp, P ee 2-K(4> pl \Hp Q \4> Q )(4> Q \H Q p\dp+) , (24) 

we thus obtain from Eq. (fT9]) 

rpr _ TpipjE) 

1 P'P — 7 -tit?\\ ' {*°) 

2tt [E-a- A(E) + 
The S operator is related to the T operator by 

S = l- 2mT . (26) 



The matrix elements of S in the (P, E) basis are 

(Xp<(E')\S\xp{E)) = 6(E - E') SPIP {E) = 5{E - E'){S P , P - 2niT P , P (E)) . (27) 

S may be expressed as the sum of background and resonant parts: S — S b + S r , and similarly s P i P — s pip + s P , p . 
The background S operator S b = 1 — 2niT b describes scattering in the open channels in the absence of coupling to 
the closed channels. The resonant part of s p > p is 

-iT P , P (E) 



S P'P - ~ 77Z7-, 7t(eT ■ ( 28 ) 

2 



E-ex-AiE)^ 1 ^- 



Using Eq. (|T4|) , the adjoint of Eq. p2|) . and Eq. (|22|) . one may show that 

= E ( 6 P'P" - ^i( XP ,\V PP \(f> P „)) (4> P „\. (29) 
p" 

The quantity in parentheses is s b pipll , hence from Eq. (|24|) we find 

T P , P = J22™ P , P 4<j>+,,\H PQ \<fi Q )(<p Q \H QP \ ( j>+) . (30) 
p" 

Since S b is unitary, 

\T PP \ = 2i:\(<j }Q \H QP \cp + p )\ 2 (31) 

and r = |Tpp|. Eq. (|28|) shows that the resonance gives rise to a peak in all open channels, with the magnitude 
p 

of the resonant scattering in channel P determined by |Tpp|. The position and width of the resonance are the same 
in all channels, and the position E r satisfies the equation E r — ei — A(E r ) — 0. If the resonance is sufficiently narrow 
and r does not vary significantly, then its width is approximately T(E r ). 

If the open channel space consists of a single partial wave, the single element s PP may be written as s P (E) — e 2lS , 
where 5(E) is the scattering phase shift. The background part of s P is s b P = e 2ldb , and 6 = Sf, + S r . Eq. (fSD|) takes 
the form 2Tr(cj) p \H P Q\(j)Q)((j)Q\HQ P \cj)+) = e 2i5b T. Using Eqs. fl24]) and ([25]), we obtain 



£7-ei-A- 1 
E-e 1 -A- l , 



s P {E)=e™»[ - r -2 r \ , (32) 



and find that the resonant phase shift is given by 



^)=^" V-t ( -W (33) 
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III. INTERMEDIATE SCATTERING STATE 

In this section, we apply the formalism of Sec.|TT]to a system with three channels: channel a, channel j3, and channel 
7. The state a is an incoming channel for the colliding particles, the state (3 is a distinct open channel, and the state 
7 is a closed channel. There is no direct coupling between a and 7 but the incoming collision channel is coupled to 7 
indirectly via a sequence of two couplings: V a p and Vp-y . In our motivating example of a collision between two atoms 



of different types 



13 



16j, a is the s-wave scattering state of the ground electronic state, f3 is the p-wave scattering 



state of the ground state, and 7 is a p-wave bound state of an excited electronic state. Figure [T^, illustrates this 
system. We denote projection operators onto the three channels by P a , Pp, P 7 , and the projections \a + ) = P a \ i S + ), 
\f3 + ) = P^| X I ,+ }, and I7) = J > 7 |\f ,+ ). We also define projected Hamiltonians H vp = P u HP pi where v, p — a, (3 or 7. 
The coupled equations for our system are (recalling that V ai — V ia = 0) 

(E-H aa )\a + ) = V a p\/3+) (34) 
{E-H p)\(3 + ) = V 0a \a + ) + V M \ 1 ) (35) 
(E-H yj )\ 7 )=V lP \f3+). (36) 

Because H is Hermitian, V a p = Vl a and Vg 7 = VlL. 

In our example of atom-atom collisions, the projection operators P v (y — a, (3,"/) are of the form \nime) {ntmi\, where 
n denotes the electronic states of the colliding atoms and £ and mi are nuclear angular momentum quantum numbers. 
The channel Hamiltonians H vl , are the nuclear radial Hamiltonians and the states |a + ), I7) correspond to 

functions of the internuclear radial coordinate R: ip+(R) = (R\a+), ipp(R) = (R\j3+), ipy{R) — (R\j)- The centrifugal 
potential £ v (ly + l)/2/ii? 2 is retained as part of the kinetic energy operator K vu . The states and |/3 + ) carry the 
+ sign because they are nonnormalizable scattering states with proper asymptotic forms as R 00. The state I7) is 
a normalizable bound state, which means that Rip~ l {R) decays exponentially as R — > 00. The off-diagonal interaction 
potentials V a p(R) = (R\V a j3\R) and V^ 7 (i?) = (i?|Va 7 |i?) are functions of R. Our goal is to find the dependence 
of scattering properties on the coupling strength, so we assume that V a p{R) = VL\V a p{R) and Vp 1 (R) — r2 2 Vg 7 (-R), 
where f2i and SI2 represent the coupling intensities and V a p{R) and Vp 7 (R) are functions independent of the coupling 
intensities. In our example of Li-Cs collision in electric fields, V a g can be approximated by a Gaussian function, 
whereas fii is proportional to the strength of the electric field [161 ]. 

The P subspace defined in Section HT1 comprises two open channels, a and (3, and the closed channel subspace Q is 



10 



the 7 channel: 

( 



Hpp — 



H aa V a f3 

\ V(3 a Hpfj 



H, 



( o N 



PQ 



(37) 



(38) 



\Vp~f J 



(39) 

H QQ = (H, n ) . (40) 

The final solution can be described in terms of Green's operators containing various parts of the Hamiltonian. The 
following section provides useful notation and identities for the Green's operators for this problem. 

A. Green's operator notation and identities 

Here the symbols v and p are used to denote any of a, f3 or 7. We denote the Green's operator for the outgoing 
wave in the isolated ^-channel by 

GHE)^^-. (41) 



L vv 



This operator satisfies the equation G+(E){E — H vv ) = 1„, where 1„ is the identity operator in the Hilbert space of 
channel v. Similarly, G~(E) represents the incoming wave Green's operator, while G l/ (E) is a Green's operator with 
neither + nor - specified. It is used to describe closed channels, where G+(E) = G~(E). 
We make use of Green's operators of the form 



G ^ {E) = E+-H vv -V vp Gt{E)V pv ' ' 42 1 



where the + superscript on the symbol Q^(p+) indicates the presence of a + superscript on the E in the denominator, 
and the + following the index p in parentheses indicates the + superscript in the V vp G+V pv term. We will use similar 
symbols where either or both of these superscripts may be replaced with a minus sign or with no superscript. The 
argument to a Green's operator is assumed to be E, unless specified otherwise. 
A useful identity is 

S^)Vu P G^ = Gi;V vp g+ v+) , (43) 
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which we derive starting from 

V up (l p - G+V pu G+V up ) = (1„ - V up G+V pu G+)V up . (44) 
Writing t p — Gp(E — H pp ) and t v = (E — H VV )G+ , we obtain after factorization 

V vp G+(E - H pp - V pu G+V vp ) = (E — H vv - V vp G+V pv )G+V vp , (45) 



which yields Eq. (|43|) after multiplying on the left by G„r p+ \ and on the right by Gp/ v+ \ ■ 
Another useful identity is 

Gt (p) =G + V + GtV vp G + p V pv g+ (p) , (46) 



which is the general linear operator identity A \ B = 4 + \B A ^ B applied to 



[P+Y 



B. The Green's operator E+ -h pp 

To apply the method of Section [II] we must evaluate the Green's operator E+ ^ H , which enters Eqs. (fl"2"|) and 
l3[) . We do this by solving the inhomogeneous Schrodinger equation 



(E - H PP )\b) = \c) , (47) 

where 



'c Q } 



(48) 



V \cp) 

is a given arbitrary vector with no components on the null space of E — Hpp, which is necessary for a solution to 
Eq. (j4"T|) to exist. Every solution to Eq. (|4T|) corresponding to the correct boundary conditions is of the form 

M = M + E+-H PP {C) > (49) 
where \b') is a solution to the homogeneous equation. Using the identity 

( K) \ 

\b)=\ , (50) 

V \h) J 

Eq. (T47|) can be rewritten as 

(E-H aa )\b a ) = V af) \bp) + \c a ) (51) 
(E-H P0 )\bp) = Vp a \b a ) + \c ). (52) 
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Applying, respectively, G+ and Gt to these equations we obtain 

K) = K) + G+V a p\b f3 )+G+\c a ) (53) 
\bp) = \b'p) + G+Vp a \b a ) + G+\cp), (54) 

where \b' a ) and \b'p) satisfy the equations (E — H aa )\b' a ) = and (E — Hpp)\Vp) — 0. Substituting these equations 
into the right-hand sides of Eqs. (I51|) and ([52]) and rearranging we obtain 



(E - H aa - V a pG+V 0a )\b a ) = V a p\b'p) + VcpG+lc'p) + \c a ) (55) 
(E - H fiP - V Pa G+V aP )\b p ) = V 0a \b' a ) + V 0a G+\c' a ) + \c p ) . (56) 

Multiplying Eqs. f|55[) and (|56p by G^ip+s and £/g( a _|\ respectively, and rewriting the system in matrix form, we obtain 
/ 



(57) 



K) = Gtw Gi { p +) V a ?G+ \c a ) + \b'L) + G + a{fj+) V af3 \b' P ) 

VIM/ \Gt [a+) v Pa G+ g+ {a+) )\\cp)J \\b^) + g+ a+) vp a \b' a ) 

where {E-H aa -V a pGjV 0a )\b'^) = and {E-Hpp - V 0a G+V a p)\b'^ = 0. This is of the form of Eq. (gg]) - the reader 
may verify that the last term on the right-hand side is a solution to the homogeneous equation - and we conclude 
that 



E* — Hpp 

which can be verified by matrix multiplication with E — Hp p 



1 _ ( ^ ( p +) (E) g+ (0+) (E)V a0 G+(E) 

g^ a+) (E)V0 a Gt(E) g; {a+) (E) 



(58) 



C. Background scattering in the open channel 

Before including the coupling to the Q-channcl, it is necessary to obtain the background scattering properties in 
the P-channcl. These solutions are written in terms of scattering solutions for the isolated a- and ft- channels. A 
free scattering state for the a-channel \ao(E)) satisfies the equation {E — K aa )\ao) — 0. Similarly, a /3-channel 
free state \fto(E)) satisfies (E — Kp )\fto) = 0. In the example of atom-atom collisions with the a and ft channels 
corresponding to s and p partial waves, (R\a {E)) — jQ^nE/h 2 ) 1 / 2 ^ and (R\ft (E)) = ji^/iE/fP^^R), where jt 
are Riccati-Bessel functions. Eq. f| 12[) for the isolated a and ft channels gives 

\<t>+) = \a ) + G+V aa \a ) (59) 
\<t>+) = \fto) + G+Vpp\ft ), (60) 
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where \<fi+(E)) is a scattering state solution for the isolated a channel, satisfying the equation (E — H aa )\<fi+(E)) = 0, 
and \(j>t(E)) is a scattering state for the isolated /3-channel and (E — H/3p)\<f>t(E)) = 0. 

For the P channel background scattering problem, we seek states \<j>p) satisfying (E ~ Hpp)\<fip+) = with 
appropriate boundary conditions. A basis for free states of the two-dimensional P channel is given by the vectors 

/ 



IXp) 



\X P ) 



l a o) 






VIA)) 



(61) 



(62) 



Equation (fT2|) relates the scattering state \4>%) to a free state \xp) in the space of the P channel. Using Eq. (|58|) for 
the Green's operator E +} Hpp , as well as Eqs. (|37|) . (|43|) . |46|) . (|59|) . and (|60|) . we solve Eq. (fl~2"|) for the free states 
(|6T|) and (62|) to obtain 



4 + ) = 



( Kg) 1 


■I 


v i^g> y 




f Kg) ^ 


■I 


v i/?b + g ) y 





(63) 



(64) 



s^ +) v aP \tf>p 

The first solution, \4>p+), represents a state with incoming flux only in the a channel, while the state \4>p + ) corresponds 
to incoming flux in the [3 channel only. 

For elastic scattering in the a-channel, the background part of the on-shcll T-matrix clement, from Eq. (fl~5|) . is 



TL = (Xp\VppW p + ) 



(a \V aa \a£ G ) + (a \V a p\f3£ G ) 



(65) 



D. Transition amplitudes 



The effective potential in the P-channel, defined by Eq. ©, takes the form 

f,r t r \ 



V eS (E) 



\ Vp a Va/5 + VfrG^fyVtf ) 
Eq. (115|) gives the resonant part of the on-shell T-matrix element for elastic scattering in the a-channel, 

1 



(66) 



T, 



.... - {4> P Wpq 



E-Hr 



-H QP \y+) 



(67) 
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The explicit form of the operator in this expression is 




so 



HP *E^ H ^= I (68) 



TL = (Pbg\ V Pi G i V i0\P + )- ( 69 ) 



Because the resonance is due to a single isolated bound state in the 7-channel we may approximate 



G,(E)^^21, (70) 



and 



Eq. (fTU)) becomes 



£-e 7 



rpr _ (/?BGl^7l^7)(^l^l/? + ) (7U 

E-e 1 ■ [ ' 



\a+) = \a+ Q ) (72) 

E — e 7 



To find the matrix element (0 7 |V r 7 ^^^ a+ j Va 7 |0 7 ), we multiply Eq. (|73|) on the left by (0 7 |T^y/3, and isolate 

(h\V^\P + ) = ( ff"fH iTT • (74) 

Substituting this result into Eq. (f7Tj) . we arrive at 

T r (/jcj^M^kfe) (75) 

E-e 1 - (^ 7 |^ 7 /3^ (Q+) V> 7 |(/) 7 ) 

In order to evaluate the matrix element {'t'-yV^-y i3Q^ a+ ^p-y\'t ) -y) appearing in the denominator of Eq. (|75p . we need to 

expand the operator Q% a+ \ in a basis of eigenstates for the operator -ff/3/3 + Vp a G^ V a p, which is unknown. However, 

since G~pi a+ \ is the /3/3-element of E+ _} Hpp , we can use an orthonormal basis of eigenstates for E — Hpp, given by 

Eqs. dMJ-lIM]), to represent 

a+ (F) r, Fl \& G (E')m G (E')\ \{3> B + G (E>))(P'+ G (E>)\ 

^(a+) W ~ J o dE ^+3^ + ■ ( ?6 ) 

Substituting this expansion into the matrix element {<j>-y\V~ l pG~p^ a+ j(E)Vp~ l \4>~ ( ) in Eq. (|75p . we obtain 



Re(</» 7 |y 7 ^+ (ct+) (^)V^ 7 |0 7 )=P 







A(E) , (77) 
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where V indicates the Cauchy principal value integral. With the use of Eq. (|22|) , we obtain 

h^{^\V^g+ [a+) {E)V^) = -7T (\{^\V lP \P+ G {E))\ 2 + \(^\V^\f3'+ G (E))\ 

_ m) 

2 ' 

Using this definition, Eq. (fT7| , and the definitions 



(78) 



T aa {E) = 2tt(P bg (E)\V^)(^\V 10 \(3+ g (E)) (79) 
T p(E) ee 2n(^ G (E)\Vp^)(^\V J p\(3'+ G (E)) , (80) 



T -^) = —7 5^ • <^ " 

2vr (£-e 7 -A(£)^^ 



we obtain 



c 7 ; -r 2 

One may confirm that the definitions in Eqs. (Tf7 |) -([80 |) agree with those in Eqs. (f2"T|) . (123| and ([24]) . As discussed in 

Secini T(i?) = |r oa ,(£?)| + \Tpp(E)\. As the coupling amplitudes Oi, f2 2 approach zero, Q„m+) ~ ¥ ^J' ^s(a+) 

From Eqs. (J5SJ) and §3§ we find that |/3+ G ) = O(fii), while from Eqs. © and ((MJ), we obtain \/3'+ G ) = 0(1). Hence 

A(J5) = o(n 2 2 ), r aa = o(n?n§), Tpp = o(n§) and r = o(n§). 

IV. INTERMEDIATE BOUND STATE 

In this section, we consider the problem in which the intermediate state is a bound state of a closed channel. The 
system is illustrated by Fig. [TJd. As in Section IIII1 we denote the scattering state solution in the full Hilbert space 
by \^ + (E)) and define projection operators P a , Pp and P 7 . In contrast to Sec. IIIIl the projection of \^ + ) onto the 
/3-channel is a bound state PpY$> + ) = \0). Our system is then described by the coupled equations, 

(E-H aa )\a+) = V a p\/3) (82) 
(E-Hpp)\f3)=Vp a \a+)+Vp 7 \ 1 ) (83) 
(E-H yj )\ 7 )=V J p\{3). (84) 

The bound states \4>p) and |<^> 7 ) in the j3- and 7- channels satisfy the equations (ep — Hpp)\f3o) = and (e 7 — -ff 77 )|7o) = 
0. We assume that these states are well-separated from other energy eigenstates in their respective channels, both 
bound and continuum. Note that ep and e 7 do not need to be well-separated from each other. 

Equations |82 p -(|84 p can be solved as described in Sec. [Til with the a-channel taken as the P-channel and the (3- and 
7-channels together forming the Q-channel, and the Q-channel may be diagonalized as described in Ref. 



m (p- 157). 
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However, it is simpler to apply the technique described in Sec.|TT]twice, first to eliminate the closed 7-channel, reducing 
the problem to the a- and f3- channels, and then to eliminate the closed [3 channel, which yields an effective equation 
for the a-channel. 

Inverting Eq. (|84p by acting on the left with G 7 and substituting the result into Eq. gives 

(E - H pp - V 07 G 7 V 70 )\P) = V Pa \a+) . (85) 

This eliminates the 7-channel from the problem and reduces the three-state problem to an effective two-channel 
problem. Applying GpM to Eq. (|55|) and substituting the result into Eq. gives 

(E - H aa - V a pQ m Vp a )\a+) = . (86) 
The effective potential experienced in the a-channel is therefore 

V cS {E) = V aa + V aP g m (E)Vf, a . (87) 
Rearranging Eq. ([55]) and multiplying on the left by gives 

\a+) = \4>t) + GtV a0 g m V Pa \a+) , (88) 

where |</>+) is, as in Sec. IIIIl a scattering state solution for the uncoupled channel. The incoming-wave stationary 
state \<j>a) is related to the free state |ao) by 

(<l>-\ = {ao\ + (ao\V aa C%. (89) 

The on-shell T-matrix element for elastic scattering in the a channel is 

T aa = (a \V cS \a + ) 

= {ao\V aa \(j>+) + {4>~ \V a pG 0M V 0a \a + ) 

= T b aa + T r aa . (90) 
Using the isolated state approximation for \4>p), we represent the Green's operator GpM by 



to obtain 



T r _ (<t>a\VaP \M (0/3 \V0 a \<X+) fq2 s 

aa E - e - (falV^G^plfa) ■ [ ' 



The matrix elements {(f>p\Vp a \a + ) can be evaluated by multiplying Eq. (j88[) on the left by (4>p\Vp a 

E—£p — {<t>p\Vp- l G~/V. r p\tf>i}) 

Substituting this result into Eq. (f§2"| yields 



E — e/3 - {^p\Vp a G%V a p\<j>p) - {(j>p\Vp 1 G 1 V 1 p\<f)p) 
Using the isolated state approximation for |^ 7 ), the Green's operator G 7 can be written as 

Gy(E) 

and hence 



E-e 1 ' 



E- 1 — E — e^y 

The states ^^{E')) form an orthonormal basis for the a-channel so we may expand as 



G+(E) = / dE> 
Jo 



E+ - E' 



to obtain 

Re (MV, a Gt(E)V a pm =V[ , A(E) . 

where V indicates the Cauchy principal value integral, and 

Im (<f>p\Vp a G+(E)V a p\<f>p) = -K\{<t>p\Vp a \<t>+{E)) | 2 = -- 



, 2 _ TO 



Inserting Eqs. ([MD, JMH and O into Eq. (JUJ) we have 



rpr 



2tt - ep - A(E) 
From Eq. (|29|) . 



2 _E-£, 



A"| = (1 - 27rf<a o |y aa |0+))^+| = s b M\ 



and with Eq. (| 1 OOf) we obtain 
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The S-matrix element for the a-channel s a is defined in Eq. (|27|) . From the discussion in the last paragraph of Sec.HIl 
it satisfies the equations s a — e 2l ( s b+s r ) _ s fc _ 2niT^ a , therefore 

e 2 * = 1 - 2me-™» T r a = e^){E - A 



and the resonant phase shift is 



S r = arg 



(104) 



Feshbach resonances correspond to poles of the S'-matrix, and therefore of T^ a , which, if A ^ 0, occur when 

(£-e T ) [E-ep-A(E) + 1 -^^- ) -A = 0. (105) 



Denoting the roots of this equation by Ex and E 2 , we see that in the limit of A — > 0, Ei — > e 7 and E 2 — tp — A{E 2 ) + 
' r ^ 2 ' ) — * 0. For A = the a and /3 channels decouple from the 7 channel and the pole at E\ is not present. The 
quantities A[E) and T(E) usually vary slowly with E so there should be two resonances, one associated with the 
bound state 1 0/3 ) and one with 7 ). If we assume that A(E) and T(E) are independent of energy, then, to first order 
in A, the poles occur at 



2A 

Ei — e 7 : fr 

e fj + A - f - e 7 

E 2 =ep+A-- + - iT ■ (106) 

2 e + A - — - e 7 



The poles move in opposite directions with increasing A. In particular, the pole at E\ moves away from the real axis, 
which results in a non-zero resonance width. In the limit of A — > 00, 

Ex,2 = e 7 + £ /3+ A ~f ± A l/2 (1 + ; (107) 

where the + root is 2?i if e 7 > ep + A and E 2 if e/3 + A > e 7 . 

The coupling between channels a and /? causes a displacement in the position of the {(ftp) resonance to + A 
when A = 0, hence degeneracy of the resonances occurs when e 7 = + A. The resonance poles are then located at 
E h2 = e 7 -ir/4±(A-(r/4) 2 ) 1 / 2 . If 4A 1 / 2 < T there are two resonances at the same energy with different widths, if 
4A 1 / 2 > r there are two resonances at different energies, both with width T/4, and when T = 4A 1 / 2 there is a single 
pole of order two. 

Fig. [2] illustrates motion of the resonance poles in the complex energy plane as A varies for a system with constant 
A and T, for a nondegenerate case. 
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FIG. 2: Trajectories of the roots of (E — Oy)(E — A — + jE) — A in the complex _E-plane as A increases. The calculations 
are performed for (eg + A)/T = 0, e 7 /r = 0.1, and a finely-spaced grid of A/Y values. Symbols indicate the points on each 
trajectory corresponding to A/T = 0, 0.1, 1 and 10. 

The left-hand side of Eq. (|103p can be factored as 

E _ £l(E) + !E^)V E _ 6(E) + 5iffiV (108 ) 



where the functions £i(E), 62(E), T\(E) and T 2 (E) are real- valued for real E and satisfy 

El(E) + 6(E) _ , ( lss+hss. ) = e , + Ec + A(E) _ m 



Fi(E)r 2 (E) : ( £ 1 (E)T 2 (E) + S 2 (E)T 1 (E)\ ( , A/1 ™ A .-^ T ( E ) 



£t(E)£ 2 (E) " ; - * { " ' zv ; 2 zv ' " ; J = e 7 ( £/3 + A(£)) - A - z^-^ . (109) 

Each pole of T^ a corresponds to a zero of one of the factors and the resonant phase shift from Eq. (|104[) can be 
expressed as the sum of phase shifts for each resonance, 

S r (E) = 51(E) + 81(E) = arg (e - S 1 (E) - + arg (e - £ 2 (E) - l I^Tj . (HO) 

As the coupling amplitudes f^i and Sl 2 vanish, A = 0(fi|), T — 0(£l\) and A = 0(Vl\). If and Vt 2 are small 
compared to the distance between the zeroth order resonance positions, ep + A — e 7 , the energy shift for the resonance 
due to \(f>-y) is Re(£i) — e 7 = 0(il 2 ), and the resonance width is Im(£i)/2 = 0(f^£7|). The resonance due to |<^g) 
has an energy shift Re(-E 2 ) - £/3 = O(Qf) and width Im(E 2 )/2 = O(fif). 
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V. MULTIPLE INTERMEDIATE BOUND STATES 

Here, we generalize the results of the previous section to the problem with several intermediate bound states, each 
in a different closed channel, coupled in a chain to the resonant bound state. We use the same notation as before, 
except that it is now convenient to label channel indices with numbers rather than letters, and to place the open 
channel last in the indexing scheme. Our system is described by the coupled equations, 

{E-H n )\l)=V 12 \2) 
{E-H 22 )\2) = V 21 \l) + V 32 \3) 



(E - H nn )\n+) = V n - hn \n - 1} , (111) 

where channel n is open and contains the incident flux and channels 1 through n—1 are closed. We assume that each 
closed channel k contains a bound state \<pk) satisfying (E — Hkk)\<t>k) = 0, and that each such state is well-separated 
from other states in its channel. We adopt the procedure of the previous section to repeatedly remove the last closed 
channel. After eliminating channels 1 and 2, we obtain the following equation for channel 3: 

(E - H 33 - V 32 g 2(1) V 23 )\3) = y 34 |4) . (112) 

For the more complicated Green's operators that appear in solving this system, such as that for the pseudo- 
Hamiltonian on the left-hand side of Eq. (|112p , we extend the notation defined in Sec. IIIIAl with a recursive definition. 
If at with k = 1 to n are channel indices, 

^ B (a n _ 1 (-(a 1 )..-))(- B ) = p — 77 T7 a — " Jew ' ( 113 ) 

The inverse of the operator in Eq. (1112[) is denoted by £3(2(1)) (E). Repeating this process we obtain for channel n an 
effective Schrodinger equation 

(E - H nn - K,n-l£ri-l(n-2(---(l)---))K-l,n)l n+ ) =°, (114) 

with the effective potential determined by 

V e ff(E) — V nn + ^, ra -l5n-l(n-2(---(l)--- ))(E)V n -l,n > (H5) 
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where the operator G n -i(n-2(---(i)---)){E) is 

Qn-l{n-2(-(l)---)){E) 



E — H n -l,n-l — V n -i,n-2 V n -2,n-l 

1 

E — H n -2,n-2 — V n -2,n-i Ki— 3,n-2 



(116) 



E — H 2 2 — V21 -5 5— V12 

ii — Til 1 



Using Eq. (|95|) repeatedly we obtain, for fe = 1, .... n — 1, 

g^-K...))^)^!^)^!^^, (117) 

where P k (E) are polynomials in _E satisfying the recursion relation 

Po(E) = 1 
P 1 (E)=E-e 1 

P k (E) = (E - e k )P k ^(E) - A k P k _ 2 (E) , fc>2, (118) 

where 

A fc = ]<0fc|Vfc,*-i|0fc_x>r- (119) 

This gives the implicit equation for |n + ), 

= + G n 14,, l -i|0„-i)<0„-i|K-i,„|n + )§^ i , (120) 
from which we can obtain the resonant contribution to the T-matrix element for elastic scattering in the open channel, 

Kn = (ClK,TO-i|^-i}(0TO-i|K,»-iIn + )^. (121) 

J n— 1 

We evaluate the matrix element (</>n-i |U„. n _i |7i + ) by multiplying Eq. (|120p on the left by {4> n -i\V n -i in : 

(0n-l\Vn-l,n\(l)i) 

1 — {4>n-l\Vn-l,nGnV n ,n-l\4>n-l) p*Zi 



-iiK,„-iin + ) = — — z . . — • ( 122 ) 



Using this result, Eq. p21[) yields 

(^-|y njn _ 1 |^_i)((^_i|y n _ 1 , n |^+>p n _2(s) 



' r "" ii: ' ' Pto-i(-E) - (^„_i|V r n -i, Tl G! B (S)V n , Tl -i|^,-i)P w _2(B) 

r(s)P„- 2 (E) (123) 



2vrQ„_ 1 (i;) ' 
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with Q n _x(E) defined by 



Qn- X {E) = P„_! 



(E) - A(E) 



iT(E) 



2 



) 



Pn-2(E), 



(124) 



where in these last two equations T(E) and A(i?) are as given by Eqs. and (|99p but with a replaced by n and (3 
by n — 1 . 

Resonances are associated with the roots of the equation Q n ^i(E) = 0, which correspond to poles of the S-matrix. 
The root which approaches efc (or e„_i + A(e„ — 1), when k = n — 1) as the coupling strengths Ai tend to zero is Ej~. 
The real part of the root gives the resonance energy, and the width is — 2Im£ , fc. A real root of Q n _i, for example 
when one of the couplings Ak = 0, must be a root of both P n ~\ and P n -2, and by Eq. (|123p does not give rise to a 
pole in T^ n . Physically, the resonance width approaches zero as the root approaches the real axis. A factorization 



In order to elucidate the properties of the resonances induced by indirect coupling, we consider the case with 
constant A and T and all Ak equal to a constant A. Q n _i{E) is then a polynomial in E and A, which may be solved 
numerically. £k ~ i^k /2 = Ek in Eqs. (| 1 25[) and (| 1 26[) implies that the roots of the polynomial sum to A — ?T /2 + £ fe ■ 
Fig. [3] illustrates the trajectories of the roots of an eighth-degree polynomial Q%(E) in the complex energy plane as 
A increases. The energies of the bound states are assumed to be closely and regularly spaced. For all roots, 
Im Ek < when A > 0, and ImEk approaches a constant as A — > oo. We have solved a range of other examples, 
varying the ordering of ek, the regularity and magnitude of their spacing, and the number of bound states n — 1, and 
found that this is a generic property. If all roots lie in the ImE < half-plane, the resonance widths sum to T and 
every resonance must have width less than T. It can be proven that the roots become real only when A = so each 
root must remain on one side of the real axis for A > 0. From Eq. (|124p . a real root of Q n -x must be a root of both 
P n -i and P n -2- From Eq. (|118p . APk-2 = [E — €k)Pk-i — Pk for k > 2, hence if A ^ 0, a common root of P n _i and 
P n -2 must be a root of all Pk, k < n — 1, including Pq. However, this would imply that Pq has a root, which is not 
possible since Pq = 1. We conclude that in a solution of Eq. (|118[) with A ^ 0, no two consecutive Pk may have a 
common root. As a consequence, there can be no real roots of Q n -\ unless A = 0. 

An analysis of the resonance widths shows that the width of the resonance at E\ due to the last bound state \4>i) 
decreases rapidly as the number of intermediate bound states increases (Fig. 2]). 



Q„_i(B) = l[E-£ k (E) + ^ 



(125) 



exists, where £k{E) and Tk(E) are real for real E, Ek — £k(Ek) ~ lTk ^ k '> an d 




(126) 
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FIG. 3: Trajectories of the roots of an eighth-degree polynomial Qs(E) in the complex B-plane as A increases. Calculations 
are performed with all Ak = A, (es + A)/r = 0, e^/r = 0.1(8 — k) for k = 1, . . . , 7, on a fine grid of A/Y values. Symbols 
indicate the points corresponding to A/T = 0, 0.1, 1 and 10. 

VI. SUMMARY 



We have presented analytical solutions of the Schrodinger equation describing multi-channel Feshbach resonances 
mediated by second- and higher-order couplings, i.e. scattering resonances induced by the interaction with a bound 
state that is not directly coupled to the initial scattering state. The results obtained represent three distinct problems: 
(i) resonant scattering induced by a sequence of one continuum-continuum and one continuum-bound couplings; (ii) 
resonant scattering induced by a sequence of one continuum-bound and one bound-bound couplings; and (iii) resonant 
scattering induced by a sequence of one continuum-bound and several bound-bound couplings. 

The first problem is a model of electric-field-induced resonances in ultracold collisions of alkali metal atoms 
Electric fields couple s-wave collision channels with p-wave scattering states. The numerical calculations of Refs. [K; 




and 16[ showed that the s-wave collision cross section may undergo a resonant variation in the presence of electric 
fields if the corresponding p-wave collision channel is coupled resonantly with a p-wave bound state. Section IIIII 
of the present paper shows that this is a general phenomenon. Eqs. (|77|) - (jFlj) provide a general description of 
such three-state Feshbach resonances and demonstrate how the width and position of the resonances in the initial 
scattering state depend on the strengths of the continuum-continuum (s-to-p) and continuum-bound (p-to-p) couplings. 
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FIG. 4: Width of the resonance due to \<j>i) for a system with n — 1 bound states calculated with all = A, (e n _i + A)/r = 0, 
e fc /r = 0.1 (n - 1 -fe) for k = 1, ...,n-2. 

Similar resonances may arise in collisions of atoms or molecules with significant magnetic or electric dipole - dipole 
interactions. The dipole - dipole interactions couple different partial waves of the collision complex and some of the 
dipolar resonances observed in ultracold collisions of Cr atoms 
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28 1 may have a three-state character. Eq. (fSTj) 



of the present work provides a general form that can be used to fit the experimental data on three-state Feshbach 
resonances involving two continuum states and the theory of Sec. |IIII may be used for a refined analysis of Feshbach 



resonances in the Bose-Einstcin condensate of Cr atoms 
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The second and third problems considered in this work are relevant for collisions of complex polyatomic molecules 
with ultracold atoms or molecules. Polyatomic molecules can potentially be cooled to ultracold temperatures by 
clastic collisions in a reservoir of ultracold atoms 
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2fJ. The experimental realization of this method may 



be complicated by naturally occurring Feshbach resonances. The energy spectrum of polyatomic molecules is quite 
dense and collisions of large molecules with ultracold atoms may lead to long-lived Feshbach resonances that would 
complicate translational energy exchange and result in sticking of atoms to molecules and the formation of clusters. 
It is therefore extremely important to understand the mechanisms of Feshbach resonances in collisions of polyatomic 
molecules with atoms. 

If the molecule is sufficiently large, the atom - molecule scattering state of interest may not be directly coupled 
to all molecular states in a collision. The atom - molecule interaction potential, however, induces couplings between 
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different states of the molecule and the entire spectrum of molecular states may be coupled to the atom - molecule 
scattering state through a sequence of one continuum - bound and several bound - bound couplings. The simplest 
example of this coupling mechanism is a collision system of a structureless atom and a diatomic molecule interacting 
through the long-range dispersion interaction. The dispersion interaction couples the ground rotational state N = 
of the molecule only with the first and second rotationally excited states N = 1 and N = 2; however, the bound 
states of the atom - molecule complex corresponding asymptotically to N > 2 may give rise to Feshbach resonances 
in collisions of ground-state molecules through a sequence ofiV>2-iV = 2 and TV = 2 - TV = couplings. 

Section IfVI presents a general analysis of three-state Feshbach resonances induced by one continuum - bound and 
one bound - bound couplings. Eq. (|102[) gives the general form of the resonant variation of the T- matrix element and 
Eqs. (|105p and (|106[) show that the two bound states give rise to two resonances. Fig. [2] illustrates the dependence 
of the resonance positions and widths on the coupling strength. Our analysis shows that the scattering amplitude 
must exhibit two resonances, even if the two bound states are degenerate. The continuum-bound coupling shifts the 
resonance energy of the first bound state, and degeneracy occurs when the shifted energy equals the second bound 
state energy. For small coupling strengths, the two resonances occur at the same position, while for large coupling 
strengths, they have the same width. 

Section [V] generalizes the results of Sec. IIVI for systems with several intermediate bound-bound couplings and 
demonstrates that Feshbach resonances may occur even if the scattering state is separated from the resonant bound 
state by a sequence of several indirectly coupled bound states. The ladder character of the couplings ensures that the 
scattering amplitude exhibits a pole near the energy of the bound state. We have shown that the energy dependence 
of the T-matrix element can be written in a general form given by Eq. (|123p . The polynomials P^ in Eq. (|123|) depend 
on the structure of the molecule and the atom - molecule interaction potentials. They can be evaluated using the 
recursive procedure described by Eqs. (|1 18[) . A numerical analysis of the polynomial roots shows that the width of the 
resonance decreases rapidly as the number of intermediate off-resonant bound states increases. Resonances induced 
by high-order couplings should be ubiquitous in collision systems involving complex molecules with multiple degrees 
of freedom and our expressions and formalism can be used for the analysis of experiments on ultracold collisions of 
atoms and molecules, chemical reaction dynamics and electron - molecule scattering. 
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